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@ Electric power system simulator. 

@ An electric power systenn simulator connprising 
first calculator which calculates a generator dynamic 
characteristic, second calculator which calculates a 
load dynamic characteristic and third calculator 
which calculates power transmission condition of an 
electric power system network, then the first calcula- 
tor cyclically obtains an amplitude and a phase of an 
internal induced voltage of the generator based on 
data on an output power and a terminal voltage of 
the generator, and the second calculator calculates a 
load equivalent admittance at predetermined time 
based on data on an A.C. voltage at a load installa- 



tion point, and the third calculator calculates a cur- 
rent flowing in each element of the electric power 
system based on respective data obtained by the 
first calculator and the second calculator and an 
electric power system network impedance to obtain 
the output power and terminal voltage of the gener- 
ator, so that a condition change of the electric power 
system including the generators is serially simulated 
by alternately executing the calculations of the first 
and the second calculators and the calculation of the 
third calculator. 
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BACKGROUND OF THE INVENTION 

Field of the Invention 

This invention relates to a simulator which sim- 
ulates, in real time, dynamic phenomena of a large 
scale electric power system, thereby analyzing 
them and further performing a function test of each 
element of the electric power system. 

Description of Related Art 

Heretofore two types of simulators are well- 
known which simulate an electric power system. 
One Is an analog type which combines many of 
apparatuses each of which simulates each element 
of the electric power system, for example, such as 
shown in Kansai Electric Power Go's catalog 
"Advanced Power System Analyzer APSA", and 
the other is a digital type which analyzes each 
element by a large general-purpose computer by 
using a calculation equation. 

However, the conventional analog type simula- 
tor can not simulate the real large scale electric 
power system and further a change of condition of 
the electric power system to be simulated is prac- 
tically impossible because it requires enormous 
work loads. Explaining more concretely, the analog 
type ^sirnulatbr" simulates the electric power system 
with a miniature model which is really provided 
with the electric power system elements such as 
generators, AVRs (Automatic Voltage Regulators), 
power generation control apparatuses like gover- 
nors, transformers, loads, SVCs (Static Var Com- 
pensator), and transmission lines, thus as a scale 
of the electric power system to be simulated be- 
comes larger, the miniature model occupies a larg- 
er space. For example, in order to simulate an 
electric power system including thirty generators, 
the above APSA occupies as large area as 700m^ 
for only a simulator, thus there is a limit to an 
occupancy area for simulation of the large scale 
electric power system. Further, when changing wir- 
ings of the electric power system to be simulated, 
it is required to change connections of the simula- 
tor and this necessitates troublesome and a large 
quantity of work loads. 

While the conventional digital type simulator 
has a few CPUs, in order to analyze differential 
equations for simulating elements such as a gener- 
ator, excitation system, PSS (Power System Stabi- 
lizer), and governing system, the simulator em- 
ploys either a batch solution which solves an equa- 
tion expressing the whole system or a partition 
solution which solves serially each equation repre- 
senting each element ("High Speed Transient Sta- 
bility Calculation Method suitable for Array Proces- 
sor", Institution of Electrical Engineers of Japan. 



treatises separate volume B, article No. 59-B36, 
pp297-304, May 1984). 

However both the solutions consume much 
processing time because the batch solution is re- 

5 quired to solve a large equation at one stroke, 
while the partition solution is required to stop all 
other processing to solve one equation during its 
solving, thus both the solutions can not cope with 
the simulation, in real time, of the real electric 

10 power system. Further, the digital type simulator 
calculates the differential equation to express the 
element of the electric power system as the case 
of the analog type simulator as mentioned above, 
but it is expected that because of necessary focus- 

75 ing calculation, an analysis of a high speed phe- 
nomenon is difficult to effect even with the decen- 
tralized computer, and further an error when mod- 
eling a real phenomenon in an equation can not be 
avoided. 

20 Now, the simulation of an electric power sys- 

tem is a task of solving alternately two sets of 
equations in time series, that is, a differential equa- 
tion that calculates dynamics of each element con- 
figuring the system such as a generator, excitation 

25 system, PSS, governing system and a rotating mo- 
tion system, and a network equation that expresses 
a relation between a bus and a transmission line. 
Fig. 1 is a flowchart showing whole of calcula- 
. . tion processes of the electric power system simula- 

30 tor shown in the above "High Speed Transient 
Stability Calculating Method suitable for Array Pro- 
cessor". 

Step S41 is a process of initialization. Step S21 
to Step S23 are processes of electric power sys- 

35 tem network calculation, Step S101 to Step SI 05 
are processes for generator dynamic characteristic 
calculation and Step S45 and Step S46 are pro- 
cesses proper to this technique. 

Procedure of this flowchart is as follows. 

40 First, after initialization in Step S41, calculation 

relating to a generator bus (Step S21) and calcula- 
tion relating to a non-generator bus (Step S22) are 
executed respectively. In this flowchart, the electric 
power system network calculation is performed 

45 considering a so-called constant power characteris- 
tic of consuming a constant power rather than a 
non-linear characteristic of a load, for example, an 
applied voltage, thus a solution is obtained by 
repeated calculations as shown in Step S21 to Step 

50 S23. That is, focusing or non-focusing is decided in 
Step S23. and in the case of non-focusing, the 
process is returned to Step S21 . 

When focusing is decided in Step S23, a gen- 
erator output power and terminal voltage are ob- 

55 tained. According to decisions* in Step S45 and 
Step S46, when focusing for the second time or 
more, the process ends, and in other case 
(focusing for the first time), a generator dynamic 
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characteristic is calculated using information on 
generator output power and terminal voltage in 
Step S101 to Step S105. That is, calculation of a 
rotating motion system (Step S101), calculation of 
PSS (Step S102), calculation of a excitation system 
(Step S103), calculation of a governing system 
(Step S104), and calculation of a generator equa- 
tion (Step S105) are respectively performed. Using 
data of an amplitude and a phase of a generator 
internal induced voltage obtained from these cal- 
culations, the process is returned to Step 821 to 
perform again the eiectric power system network 
calculation. 

In addition, in the flowchart of Fig. 1, the pro- 
cesses from "START" to "END" are perfornned in 
one time step of integration, and by cyclically re- 
peating the process of every one time step of 
integration, a change of condition of the electric 
power system is serially calculated. 

Now, in the calculation of the rotating motion 
system in Step S101, a generator rotor position 0 
Is calculated by introducing a phase 5 as a variable 
which meets an equation 

0 = 400 t + 5 

as shown, for example. In "Analytical Theory of 
Electric Power System" (pp 294-299. written by 
Sekine, published by Denki-Sholn. in 1971), where 
ct) 0 is a rated angular velocity In an A.C. system. In 
a usual A.C. system, all generators are rotating at 
an almost rated angular velocity, thus it is known 
that the phase 5 remains within the range of - 
360n* to +360n° and there is no problem. How- 
ever, when one generator causes, for example, 
step-out to increase a rotational speed or when the 
generator is in stopped condition, the phase 5 
rapidly increases or decreases. In this case, there 
arises no problem In pure theory, but a problem 
arises in a digital numerical calculation. 

In a computer system, a numerical value is 
generally expressed with "A* 10^", whereby the 
system usually handles numbers each of which 
combining a mantissa A and an exponent B. But 
hexadecimal numbers are actually used in the 
computer system, however, for simplification, the 
numerical calculation will be explained with decimal 
numbers. 

For example, in the computer system handling 
a mantissa A of four digits and an exponent B of 
two digits, a decimal number "58325" is expressed 
as 5833"10'' and this is handled as a number 
"583301". Therefore, either the case where the 
phase 5 is equal to "58325" or equal to "58333", 
both the numbers are handled as "583301 " in the 
computer system. Thus, the difference between 
both the numbers is to be "8" in pure theory, but 
becomes "0" in the computer system. Therefore, 
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when step-out or stop of a generator lasts long, the 
conventional digital type electric power system 
simulator sometimes handles data beyond the 
range of significant figures of the computer system 
6 as mentioned above and becomes unable to per- 
form precise calculations and can not simulate the 
system - 

As mentioned above, when simulation of con- 
dition of long-lasted generator step-out, simulation 
10 of a system including a stopped generator or sim- 
ulation of a system including two different fre- 
quencies A.C. systems connected with a D.C. 
transmission system are performed by the conven- 
tional digital type electric power system simulator, 

75 the value of the phase 8 becomes unusually large 
and falls within the range of a numerical calculation 
error of the computer system and this sometimes 
disables the computer system from calculation, that 
is. the simulation. 

20 Further, in the conventional electric power sys- 

tem simulator, as shown in the flowchart of Fig. 1, 
repeated calculations are employed for the electric 
power system network calculation, thus extending 
an operation time, and a number of times of repeti- 

25 tion until focusing differs depending on other con- 
dition, thereby causing a problem of an inconstant 
operation time. Therefore, it was difficult to per- 
form, for example, the so-called real time simula- 

tion which advances calculation in a constant, time 

30 interval. Further, in the so-called threshold running 
like the case of a voltage instability state, there are 
problems such that focusing ability is degraded to 
increase the number of times of calculation repeti- 
tion and sometimes the case of non-focusing ar- 

35 ises. 

Further, in the conventional electric power sys- 
tem simulator performing the procedure of the 
flowchart of Fig. 1, all processings are executed by 
one calculating unit, thus time management, that is, 

40 processings in each predetermined time interval 
are relatively easily performed, but a time neces- 
sary for calculation is extended. As processings of 
Step S101 to Step Si 05 are required to be particu- 
lariy performed on each generator, as a number of 

45 generators increases, a time of processings of Step 
S101 to Step S105 is^ more extended. In addition, 
there are problems such that, as the number of 
generators increases, the time management be- 
comes more difficult and the real time simulation is 

50 disabled. 

Now, as an application example of the afore- 
mentioned decentralized type computer, the inven- 
tion disclosed in Japanese Patent Application No. 
1-245493 (1989) is proposed as a digital type sim- 

55 ulator capable of high speed processing. . 

In the invention disclosed in the Japanese Pat- 
ent Application No. 1-245493 (1989), a decentral- 
ized type computer having a plurality of CPUs 61 

3 
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is used as shown in Fig. 2 and a different equation 
expressing each element of the electric power sys- 
tem such as a generator, excitation system, PSS or 
a governing system is assigned to each one of the 
plurality of CPUs 61. Then each calculation proces- 
sings divided into each element are decentralized- 
processed by each CPU 61, thereby the load of 
each CPU is alleviated, thus the calculations can 
be processed at a high speed as a whole. 

By the way reference numeral 62 denotes a 
communication line connecting the space between 
respective CPUs 61 . 

Fig. 3 is a flowchart showing procedures relat- 
ing to one bus which is one of the elements of the 
electric power system in the simulator of the inven- 
tion disclosed in the above Japanese Patent Ap- 
plication No, 1-245493 (1989). By the way the bus 
includes a generator bus which is directly con- 
nected to the generator and a n on -generator bus 
which is not directly connected to the generator. 

As shown in Rg. 3. in a task of obtaining a 
solution of the different equation performed by 
each CPU 61, a task of calculation relating to 
elecitic systems of each generator bus 31 (PSS 33, 
excitation system 34, generator 35) and a task of 
calculation relating to mechanical systems of each 
generator bus 31 (rotating motion system 36, gov- 
erning system 37) can be parallel-executed. How- 
ever the task of calculation relating to each element 
of. the electric system and that of the mechanical 
system are performed by CPUs 61 independently, 
thus data interchange is necessary between both 
the tasks. By the way processing on the non- 
generator bus 32 does not requires communication 
with a calculation task of dynamics of the generator 
bus 31. 

In this way. a network equation 38 representa- 
tive of a relation between the buses is calculated 
using results of the decentralized-processings of 
respective CPUs 61 and transferring data 39 be- 
tween the buses until focusing at each CPU 61 . 

Now. in the conventional simulator using the 
above decentralized type computer, each different 
equation expressing each element of the electric 
power system is assigned to each CPU 61, where- 
by respective calculations on each element of the 
electric power system are decentralized to reduce 
a processing time. However, on the other hand, 
generally a solution of the network equation is 
obtained by solving system condition expressed by 
a linear equation with the LU analysis method or 
Newton-Raphson method. For example, generally 
regarding a generator and a load as equivalent 
current sources and expressing the system con- 
dition with the following linear equation system, 

1 = YV 



where 

I: current 
Y: admittance 
V: voltage 

5 then the system condition is solved with the LU 
analysis method regarding the current I as a known 
quantity and the voltage V as an unknown quantity. 

Besides other some algorithms are being con- 
sidered which solve the equation with the decen- 

10 traiized type computer, but it is not right to say that 
high speed analysis capability of the decentralized 
system is best used even when any method is 
used. 

In addition, the conventional digital type sim- 
15 ulator has a problem in man-machine device. For 
example, the simulator can not output an intermedi- 
ate result of an analysis calculation of the electric 
power system and has only a function of outputting 
the analysis result in a format of a table or graph 
20 only after the end of the calculation. This originates 
directly from logic unable to stop the analysis cal- 
culation or to intermediately output the calculation 
result, but besides that it is supposed that because 
of a low speed of the analysis calculation, an 
25 operator could not feel a real time sense and he 
did not also sense necessity of the man-machine 
device except output of the analysis result 

__SUMiyiARY OF THEJNVENTION . _ 

30 

This invention was attained in view of the 
above problems belonging to the conventional 
electric power system simulator. 

The first object of the present invention is to 

35 provide an electric power system simulator which* 
can eliminate the problems originated from repeat- 
ed calculations for the system network analysis and 
can cope with the real time simulation. 

The second object of the invention is to pro- 

40 vide an electric power system simulator which can 
perform calculation for the system network analysis 
in a short time, and can execute the real time 
simulation by performing relatively correct time 
management for the calculation even when a num- 

45 ber of generators in the electric power system to 
be simulated is increased. 

The third object of the invention is to provide 
an electric power system simulator which prevents 
a generator phase from exceeding a predetermined 

50 value and also prevent the phase value from falling 
within the range of a numerical calculation error 
when the generator causes step-out. whereby the 
simulator can cope with simulation of condition 
lasting long after the generator step-out, simulation 

55 of a system including a stopped generator or sim- 
ulation of a system having two different frequencies 
A.C. systems connected with a D.C. transmission 
system. 



BNSDOCID: <EP 044571 3A2_L> 



7 EP 0 



The forth object of the invention is to provide a 
real tinne electric power system simulator which 
can quickly process focusing calculation in the 
electric power system simulator using a decentral- 
ized type computer and can simulate a phenom- 
enon of a large scale electric power system at a 
speed near the real phenomenon proceeding 
speed. 

The fifth object of the invention is to provide an 
electric power system simulator which can simulate 
phenomena of a large scale electric power system 
with easy setting and can easily changing the elec- 
tric power system and further can analyze phenom- 
ena of real machines and apparatuses which are 
difficult to analyze by a computer by making the 
best use of each advantage of an analog type 
simulator and a digital type real time simulator. 

The sixth object of the invention is to provide 
an electric power system simulator which can out- 
put not only a final result but also intermediate 
condition in real time and further comprises a man- 
machine device which is able to change condition. 

The first invention comprises first calculating 
means which calculates a generator dynamic char- 
acteristic, second calculating means which calcu- 
lates a load dynamic characteristic and third cal- 
culating means which calculates power transmis- 
sion condition of an electric power system network. 
The first calculating means . receives data - on a 
generator output and a terminal voltage to obtain 
an amplitude and a phase of the internal induced 
voltage of the generator after every one time step 
of integration, and the second calculating means 
receives data on an A.C. voltage at a load installa- 
tion point to calculate a load equivalent admittance 
after every one time step of integration, and the 
third calculating means receives each data ob- 
tained by the first calculating means and the sec- 
ond calculating means, and an electric power sys- 
tem network impedance to calculate an A.C. cur- 
rent flowing In each element of the electric power 
system by an algorithm for the system network 
calculation to obtain output of the generator and 
terminal voltage, and then a condition change of 
the whole electric power system including gener- 
ators is serially simulated by alternately performing 
the calculations executed by the first and the sec- 
ond calculating means and the calculation executed 
by the third calculating means. 

In the first invention mentioned above, a non- 
linear characteristic of the load is considered in the 
dynamic characteristic of the load, and in the case 
of the electric power system calculation, the net- 
work calculation is performed using an amplitude 
and a phase of the internal induced voltage of the 
generator, a load admittance and an impedance of 
the power transmission system, whereby the in- 
duced voltage is processed as an independent 
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voltage source and a known quantity, thus voltage 
and current values at any point are obtained from 
the network calculation, that is, a combination of 
matrix calculations, therefore the conventional re- 
5 peated calculations become unnecessary, as a re- 
sult a quantity of calculation becomes constant and 
this enables the real-time simulation. 

The second invention comprises first calculat- 
ing means which calculates a generator dynamic 

10 characteristic, second calculating means which cal- 
culates power transmission condition of an electric 
power system and communicating means which 
transfers information between the first and the sec- 
ond calculating means. Then the first calculating 

15 means receives data on an output power and a 
terminal voltage of a generator to calculate data on 
an internal induced voltage of the generator after 
every one time step of integration to output the 
data to the second calculating means via the com- 

20 municating means, and the second calculating 
means receives the data outputted from the first 
calculating means and the various constants of the 
electric power system network to calculate power 
transmission condition of the electric power system 

25 network at a predetermined time and outputs data 
on the output and a terminal voltage of the gener- 
ator to the first calculating means via the commu- 
nicating means, and then a dynamic change of the 
electric power system including generators is seri- 

30 ally simulated by performing alternately calculation 
processings by the first and the second calculating 
means. By the way the first calculating means 
includes a calculating unit which calculates the 
power transmission condition of the electric power 

35 system network and a time management unit which 
directs the calculating unit to start at every one 
time step of integration. 

In the second invention mentioned above, 
when a number of generators in the electric power 

40 system to be simulated is increased and thus a 
quantity of calculation Is increased, a new calculat- 
ing unit is added to take on calculation on one or a 
plurality of generators, whereby extending of an 
operation time can be avoided. Further, in the 

45 second invention, time management of start of cal- 
culation is performed at every one time step of 
integration, thus the time management for the sim- 
ulator becomes possible, as a result the real-time 
simulator is realized. 

50 The third invention comprises first calculating 

means which calculates a generator dynamic char- 
acteristic and second calculating means which cal- 
culates power transmission condition of the electric 
power system network. Then the first calculating 

55 means receives data on an output power and a 
terminal voltage of the generator outputted from the 
second calculating means to calculate an amplitude 
and a phase of an internal induced voltage of the 

5 
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generator after every one time step of integration to 
output these data to the second calculating means. 
The second calculating means receives the am- 
plitude and phase of the internal induced voltage of 
the generator outputted from the first calculating 
means and various constants of the electric power 
system network to calculate power transmission 
condition of the electric power system network at a 
predetermined time to output data on the output 
power and a terminal voltage of the generator to 
the first calculating means. Further when a phase 
of the first calculating means exceeds a predeter- 
mined set value obtained by multiplying 360* by 
an integer, a value subtracted from this phase 
value by the predetermined value is determined as i 
the phase, and then a dynamic change of the 
electric power system including generators is seri- 
ally simulated by performing alternately calculation 
processings by the first calculating means and the 
second calculating means. 2 

In the third invention mentioned above, after 
the phase 5 of the internal induced voltage of the 
generator was calculated as usual, a check is made 
to know Whether the phase value exceeds the 
predetermined set value or not, and when exceeds, 2. 
the phase 5 is so controlled that it remains within 
the range of the set values by adding or subtrac- 
ting the set value to or from the phase 5 , as a 
result troubles with calculation do not occur in the . 
computes system. 3^ 

In the fourth invention, when simulating an 
electric power system with a decentralized type 
computer, a simulator performs a network calcula- 
tion of one bus only on the basis of information of 
buses directly connected to the bus. Thereby the 36 
network calculation focuses fast, thus the analysis 
in real time becomes possible. 

In the fifth invention, an analog type electric 
power system simulator and a real-time digital type 
electric power system simulator are connected with 4o 
an A/D-D/A converter to transfer data between both 
the simulators, whereby large scale real-time sim- 
ulation can be performed by making the best use 
of advantages of both the analog type electric 
power system simulator and the real-time digital 45 
type electric power system simulator. 

In the sixth invention, an electric power system 
simulator of the present invention which can pro- 
cess calculations of a phenomenon of the electric 
power system at a speed near the proceeding so 
speed of the real phenomenon and outputs electric 
power system condition every one time step of 
integration during processing and comprises a 
man-machine device which can display the output 
data and perform interruption of a condition change 55 
at any time. Whereby the electric power system 
condition is displayed on the man-machine device 
in real time in a sense near feeling of the real 



phenomenon of the electric power system, and 
further time management of operation time be- 
comes possible, as a result data display and pro- 
cessing are performed as if real troubleshooting 
were executed. 

The above and further objects and features of 
the invention will more fully be apparent from the 
following detailed description with accompanying 
drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a flowchart showing the whole of pro- 
cessing steps of the conventional electric power 
system simulator, 

Fig. 2 is a schematic diagram showing a struc- 
ture of an electric power system simulator using 
a conventional decentralized type computer hav- 
ing a plurality of CPUs and of the fourth inven- 
tion of the electric power system simulator of 
the present invention. 

Fig. 3 is a flowchart showing processing proce- 
dures relating to one bus constituting a element 
of the electric power system in the electric pow- 
er system simulator shown in Fig. 2, 
Fig. 4 is a block diagram showing one example 
of a configuration of the first invention of the 
electric power system simulator. 
Fig, 5 is a circuit diagram showing a concept of 
an electric power system. 

Fig, 6 is a block diagram showing transfer func- 
tions indicating a load dynamic characteristic of 
an electric power system, 

Fig. 7 is a block diagram showing a configura- 
tion of another embodiment of the first invention. 
Fig. 8 is a block diagram showing one example 
of a configuration of the second invention of the 
electric power system simulator. 
Fig. 9 is a block diagram showing a configura- 
tion of another embodiment of the second inven- 
tion. 

Fig. 10 is a block diagram showing a configura- 
tion of one embodiment of the third invention of 
the electric power system simulator and flows of 
information at the time of calculating each ele- 
ment. 

Fig. 11 is a block diagram showing a configL u- 
tion of another embodiment of the third inven- 
tion. 

Fig. 12 (a), (b) and (c) are waveform charts 
showing changes in time series of phases of 
internal induced voltages of the generator. 
Fig. 13 is a block diagram showing a structure 
of the fifth invention of the electric power sys- 
tem simulator. 

Fig. 14 is a schematic diagram showing a struc- 
ture of the sixth invention of the electric power 
system simulator, and 

Fig. 15 is a functional block diagram of the 
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electric power system simulator shown in Fig, 
14. 

DESCRIPTION OF THE PREFERRED EMBODI- 
MENTS 

The present invention will be described refer- 
ring to tlie drawings of the embodiments. 

Fig. 4 is a block diagram showing one example 
of a configuration of the first invention of the elec- 
tric power system simulator. 

The first invention comprises generator dynam- 
ic characteristic calculator 1 as first calculating 
means, power transmission condition calculator 2 
of electric power system network as third calculat- 
ing means, and load dynamic characteristic cal- 
culator 3 as second calculating means. 

Fig. 5 is a circuit diagram showing a concept of 
an electric power system, and Fig. 6 is a block 
diagram showing transfer functions indicating a 
load dynamic characteristic of an electric power 
system. 

In Fig. 5. reference numeral 4 denotes a volt- 
age source (Voe') for an internal induced voltage of 
the generator, numeral 5, an equivalent resistance 
(conductance) G of a load, and numeral 6, an 
electric power system network. 

In Fig. 6, reference numeral 7 denotes a pri- 
mary delay system 7 for an input Va, numeral 8. a 
function which outputs "l/X^" to an input "X". and 
numeral 9, a block which multiplies any input by a 
constant Go- Therefore, in the block diagram of Fig. 
6, the conductance G is obtained concerning the 
input Va. For simplification of explanation, only the 
conductance G is treated but susceptance B can 
be also treated in the same way. 

The generator dynamic characteristic calculator 
1 performs calculation searching for a difference 
between movement of a generator rotor, that is, 
output energy of a turbine and electric energy of 
the generator and calculation for computing an 
outputted voltage, that is, a current flowing through 
a field circuit. By the generator dynamic char- 
acteristic calculator 1, an amplitude Vo and a phase 
6 of the internal induced voltage of the generator 
are calculated, and G, B of the load admittance are 
calculated by the load dynamic characteristic cal- 
culator 3. 

As can be seen from Fig. 5, all variables nec- 
essary for the electric power system network cal- 
culations, that is, an amplitude Vo and phase 8 of 
the voltage source 4, an admittance (conductance 
G and susceptance B) of the load 5, and an imped- 
ance as system information are determined by the 
generator dynamic characteristic calculator 1 and 
the load dynamic characteristic calculator 3. 

When these values are determined, voltages at 
all nodes and currents in all branches are obtained 
as is shown in, for example, "Graph Theory Cir- 



cuit" (published by Shokodo in 1974. pp109-115), 
thus by using these values, the power transmission 
condition calculator 2 of an electric power system 
network calculates an A.C. current flowing through 
5 each element of the electric power system by the 
system network calculation and can easily calculate 
an output power Pq, terminal voltage Vt of the 
generator and an A.C. voltage Vl at a load installa- 
tion point. 

10 Further, a non-linear characteristic of the load 5 

can be simulated with a dynamic characteristic of 
the load 5 as shown in Fig. 6. In a primary delay 
system 7, character T denotes a time constant and 
S, differential. Thus, assuming that Vl is constant in 

75 Fig. 6, conductance G can be expressed by the 
following equation (1). 

G = Go / Vl2 (1) 

20 Thus power P consumed in the load 5 is expressed 
by the following equation (2). 

p 

= Vl2G 
25 = Go (2) 

That is, regardless of Vo value, the consumed 
power -P becomes constant value Go, thus the 

30 constant power characteristic is simulated. 

Thus, when altemately repeating processing by 
the generator dynamic characteristic calculator 1 
and the power transmission condition calculator 2 
of an electric power system network and process- 

35 ing by the load dynamic characteristic calculator 3, 
the processing proceeds one by one time step of 
integration, and a condition change of the electric 
power system including generators can be serially 
simulated. 

40 In the above embodiment, the conventional re- 

peated calculations become unnecessary which are 
required when the power transmission condition of 
the electric power system network is calculated, 
thus the operation time becomes constant and an 

45 electric power system simulator with high speed 
operation capability can be obtained. Further, be- 
cause the operation time becomes constant, this 
simulator can cope with the real time simulation. 
Fig. 7 is a block diagram showing a configura- 

50 tion of another embodiment of the first invention. 

In this embodiment, in order to apply the elec- 
tric power system simulator of the first invention to 
an electric power system comprising three gener- 
ators and two loads, three of the generator dy- 

55 namic characteristic calculators 1 as the first cal- 
culating means 11, 12 and 13 are provided coping 
with three generators #1. #2, #3. and two of the 
load dynamic characteristic calculators 3 as the 
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second calculating means 31, 32 are also provided 
coping with two loads. Further, the power transnnis- 
sion condition calculator 2 of an electric power 
system network as the third calculating means is 
provided only one. 

Then respective amplitudes Voi, vo2 and V03 
and phases 61 , 62 and 53 of the internal induced 
voltages of the generator of the three generators 
are calculated by the generator dynamic character- 
istic calculators 11, 12, 13 and respective admit- 
tances Gi, G2, Bi and B2 of the two load are 
calculated by the load dynamic characteristic cal- 
culators 31, 32. By using these values, respective 
output powers Pqi, Pg2» Pg3 and terminal voltages 
Vn, V12. Vi3 of the three generators, and respective 
voltages Vli, Vu at two load installation points can 
be calculated with ease and at a high speed by the 
power transmission condition calculator 2 of an 
electric power system network. 

Processing functions of the generator dynamic 
characteristic calculators 11, 12. 13 and the load 
dynamic characteristic calculators 31, 32 are in- 
dependent one another and not restricted by order 
of calculation 5peration, therefore processing in 
any order of calculation is acceptable. Further, 
when a plurality of calculating units are available, 
parallel calculation can be performed by executing 
each processing function by means . of each cal- 
culating unit. Such parallel calculation can reduce 
an operation time as a whole. , . 

As mentioned above, according to the first in- 
vention, the simulator comprises the generator dy- 
namic characteristic calculator 1 as the first cal- 
culating means which calculates the generator dy- 
namic characteristic, the load dynamic characteris- 
tic calculator 3 as the second calculating means 
which calculates the load dynamic characteristic, 
and the electric power system network calculator 2 
as the third calculating means which calculates the 
power transmission condition of an electric power 
system network. The generator dynamic character- 
istic calculator 1 receives an output power and a 
terminal voltage of the generator to obtain an am- 
plitude and a phase of the internal induced voltage 
of the generator after every one time step of in- 
tegration, and the load dynamic characteristic cal- 
culator 3 receives an A.C. voltage at a load installa- 
tion point to calculate a load equivalent admittance 
at every one time step of integration, and the 
power transmission condition calculator 2 receives 
each data obtained from the generator dynamic 
characteristic calculator 2 and the load dynamic 
characteristic calculator 3 and an electric power 
system network impedance to calculate an A.C. 
current flowing in each element of the electric 
power system according to an algorithm for the 
system network calculation to obtain the output 
power and terminal voltage of the generator, and 



0445713A2 I > 



then a condition change of the electric power sys- 
tem including generators is serially simulated by 
alternately performing the calculations executed by 
the generator dynamic characteristic calculator 1 

5 and the load dynamic characteristic calculator 3 
and the calculation executed by the power trans- 
mission condition calculator 2. Therefore, repeated 
calculations necessary for the conventional simula- 
tor becomes unnecessary and a quantity of cal- 

70 culation becomes constant, as a result the real time 
simulation becomes possible. 

Fig. 8 is a block diagram showing one example 
of a configuration of the second invention of the 
electric power system simulator. 

75 The simulator of the second invention com- 

prises generator dynamic characteristic calculator 1 
as first calculating means, power transmission con- 
dition calculator 2 of an electric power system 
network as second calculating means, and a com- 

20 munication line 40 as communicating nneans which 
connects both the calculators 1 and 2. 

The generator dynamic characteristic calculator 
1 as the first calculating means performs calcula- 
tion for obtaining a difference between movement 

25 of a generator rotor, that is, output energy of a 
turbine and electric energy of the generator and 
calculation for obtaining an outputted voltage, that 
is, a current flowing-in a field circuit. -The-generator- 
dynamic characteristic calculator 1 receives an out- 

30 put power Pq and terminal voltage Vt of the gener- 
ator to output an amplitude Vo and phase S as an 
internal voltage information of the generator after 
every one time step of integration. 

The power transmission condition calculator 2 

35 of an electric power system network as the second 
calculating means receives the amplitude Vo and 
the phase 5 as the internal voltage information of 
the generator calculated by the generator dynamic 
characteristic calculator 1 and various constants of 

40 the electric power system network are used to 
perform A.C. network calculation to output the out- 
put power Pq and the terminal voltage Vi of the 
generator. 

In addition, the various constants of the electric 
45 power system network include, for example, are 
A.C. impedance of a transmission line of the power 
system network, equivalent A.C. impedance of a 
load, and winding ratio of a transformer. 

The communication line 40 transfers, between 
50 the generator dynamic characteristic calculator 1 
and the power transmission condition calculator 2 
of an electric power system network, information 
such as an amplitude Vo. phase 5, generator output 
power p G. terminal voltage V , outputted from the 
55 calculator 1 and calculator 2 respectively. 

When the processing of the generator dynamic 
characteristic calculator 1 and the processing of the 
power transmission condition calculator 2 of an 
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electric power system network are alternately re- 
peated communicating information each other via 
the communication line 40, a condition change of 
the electric power system network can be serially 
simulated and calculation equivalent to that of the 
conventional simulation can be executed by a plu- 
rality of calculators. 

In the above embodiment, system analysis cal- 
culations are divided into two parts, that is, gener- 
ator dynamic characteristic calculating part and 
power transmission condition calculating part of an 
electric power system network, thus it is possible 
to select the respective optimum systems as the 
respective calculator, as a result calculation pro- 
cessing is accelerated. Further, when a quantity of 
calculation is increased due to plural generators 
and thus new calculators are added and calculation 
for one generator or plural generators are assigned 
to one new calculators or a plurality of them, ex- 
tending of an operation time due to increase in a 
number of generators can be avoided. 

Fig. 9 is a block diagram showing another 
embodiment of a configuration of the second inven- 
tion. 

In this embodiment, in order to apply the elec- 
tric power system simulator of the second invention 
to an electric power system including three gener- 
ators #1. #2, #3. the simulator is provided with two 
generator dynamic characteristic calculators 1 of 
11, 12- By the way, a power transmission condition 
calculator 2 of an electric power system network Is 
provided only one, but includes a power transmis- 
sion condition calculating unit 2a of an electric 
power system network which Is an original function 
and a time management unit 2b. 

Further, numeral 41 denotes communication 
line connecting the power transmission condition 
calculating unit 2a of an electric power system 
network with the generator dynamic characteristic 
calculator 11, numeral 42 denotes communication 
line the power transmission condition calculating 
unit 2a with the generator dynamic characteristic 
calculator 12. 

In this embodiment, the generator dynamic 
characteristic calculator 11 performs calculation 
n .viceming the two generators #1. #2, and the 
generator dynamic characteristic calculator 12 per- 
forms calculation concerning one generator #3. 

It is assumed that the generator dynamic char- 
acteristic calculator 11 Is the so-called data-driven 
type computer which starts calculation on receiving 
output powers Pqi and Pg2, terminal voltages Vn, 
V,2 of the generators #1 and #2, respectively. Fur- 
ther, it is assumed that the generator dynamic 
characteristic calculator 12 is also the so-called 
data-driven type computer which starts calculation 
on receiving an output power Pqs, terminal voltage 
Vt3 of the generator #3. 
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Further, it is assumed that the power transmis- 
sion condition calculating unit 2a of an electric 
power system network is originally the data-driven 
type computer which starts calculation on receiving 

5 internal induced voltage amplitudes Voi. V02, V03 
and phases 61, 82, ^3 of the generators #1, #2 and 
#3. Start time of calculation of the power transmis- 
sion condition calculating unit 2a of an electric 
power system network is controlled by the time 

70 management unit 2b at every one time step of 
integration. 

Therefore, the power transmission condition 
calculating unit 2a of an electric power system 
network starts calculation when it finishes receiving 

15 of the above Information and gets a start enable 
signal for calculation from the time management 
unit 2b. The generator dynamic characteristic cal- 
culators 11 and 12 automatically start calculation 
when they output via the communication lines 41 

20 and 42 the output powers Pqi, Pg2. Pg3 and termi- 
nal voltages Vn, Via, of the generators #1, #2. 
#3 which have been calculated by the power trans- 
mission condition calculating unit 2a of the electric 
power system network. Therefore, when the time 

25 management unit 2b cyclically gives the calculation 
start enable signal to the power transmission con- 
dition calculating unit 2a of the electric power sys- 
tem network in the same time interval, the whole of 
the calculation is executed in the same time-inter- 

30 val. 

Further, when making the time interval set by 
the time management unit 2b coincide with every 
one time step of integration, variations of many 
variables during calculation processing vary in syn- 

35 chronization with real time transition, thus the so- 
called real time simulation is realized.' 

Or. the time interval set by the above time 
management unit 2b can be set in any length, thus 
the simulation at a desired speed can be realized. 

40 As mentioned above, the simulator of the sec- 

ond invention comprises the generator dynamic 
characteristic calculator 1 as the first calculating 
means which calculates the generator characteris- 
tic, the electric power system network calculator 2 

45 as the second calculating means which calculates 
power transmission condition of an electric power 
system network and the communication line 40 as 
the communicating means which communicates in- 
formation between the generator dynamic char- 

50 acteristic calculator 1 and the electric power sys- 
tem network calculator. Then the generator dy- 
namic characteristic calculator 1 receives informa- 
tion on output power and terminal voltage of the 
generator outputted from the electric power system 

55 network calculator 2 to calculate information on 
internal induced voltage of the generator after ev- 
ery one time step of integration and outputs this 
information to the electric power system network 
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calculator 2 via the communication line 40. while 
the electric power system network calculator 2 cal- 
culates the power transmission condition of the 
electric power system network at a predetermined 
time by using information outputted from the gen- 
erator dynamic characteristic calculator 2 and var- 
ious constants of the electric power system net- 
work to output information on output power and 
terminal voltage of the generator to the generator 
dynamic characteristic calculator 1 via the commu- 
nication line 40, thus by alternately performing cal- 
culation by the generator dynamic characteristic 
calculator 1 and the electric power system network 
calculator, dynamic variations of the electric power 
system Including generators are serially simulated. 
Therefore, a plurality of the calculator can be used 
depending on contents of calculation, thus when a 
quantity of calculation is increased due to increase 
in a number of generators in the electric power 
system to be simulated, calculation concerning one 
or plural generators are decentralized, thereby ex- 
tending of the operation time can be avoided. 

Further, in the second invention, the electric 
power system network calculating means 2 com- 
prises the calculating unit which calculates the 
power transmission condition of the electric power 
system network and the time management unit 
which directs the calculating unit to start. calculation 
at every one time step of integration. Thus time 
management for calculation execution can be per-- 
formed in every one time step of integration, and 
time management for the whole system becomes 
possible, as a result the real time simulation is 
realized. 

Fig. 10 is a block diagram showing a configura- 
tion of one embodiment of the third invention and a 
flow of information at the time of calculation execu- 
tion of each element 

The simulator of the third invention comprises 
a generator dynamic characteristic calculator 1 as 
first calculating means and an electric power sys- 
tem network calculator 2 as second calculating 
means. Further, the generator dynamic characteris- 
tic calculator 1, includes generator rotating motion 
system calculator 11 as phase calculating means, 
correcting means for correcting a phase S of an 
internal induced voltage of the generator 
(hereinafter to be called phase corrector 12) and 
generator excitation system calculator 13 as am- 
plitude calculating means. 

Now, an amplitude Vo and phase 5 of the 
internal induced voltage of the generator will be 
calculated by the generator dynamic characteristic 
calculator 1 . 

The phase S is calculated by, for example, the 
following equation (3). 

5 = yw o[{ (P M - P g) / M dt} - 1-0] dt (3) 
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In the equation (3), Pq is the value of an output 
power of the generator which is inputted to the 
5 generator dynamic characteristic calculator 1, and 
Wo, M, Pm are constants. 

The phase S is once calculated by the gener- 
ator rotating motion system calculator 11 and 
thereafter corrected by the phase corrector 12. 
TO That is, when it is decided in Step 12a that the 
phase 5 exceeds a predetermined value SBOn* , 
the phase 6 is determined to be "5 - 360n" at Step 
12b. On the contrary, when it is decided in Step 
12c that the phase 5 is less then -360* , the phase 
75 8 is determined to be "5 + 360n" in Step 12d. By 
performing these processings, the phase 5 is kept 
to be within the range of -360° through +360° . 

On the other hand, the amplitude Vo of the 
induced voltage of the generator is calculated by 
20 the generator excitation system calculator 13 to be 
outputted. 

The amplitude Vo and phase 5 of the internal 
induced voltage of the generator obtained by the 
generator dynamic characteristic calculator 1 are 

25 inputted to tiie electric power system network cal- 
culator 2. The electric power system network cal- 
culator 2 performs A.C. network calculation by us- 
ing information such as-an A.C- impedance of an 
electric power system network transmission line, an 

30 equivalent A.C. impedance of a load, and a winding 
ratio of a transformer to output the output power Pq 
and terminal voltage Vt of the generator. Further, 
the output power Pq and terminal voltage Vt of the 
generator are respectively inputted to the generator 

35 rotating motion system calculator 1 1 and the gen- 
erator excitation system calculator 13 of the gener- 
ator dynamic characteristic calculator 1 to calculate 
the phase 5 and amplitude Vo after every one time 
step of integration. 

40 Now, since the phase 5 is always used as an 

argument of a trigonometric function, the following 
expressions are established from the nature of the 
trigonometric function, 

45 sins = sin(5 - 360n) 
cosS = cos(5 - 360n) 

thus even subtracting or adding 360n from or to the 
phase 5 in Step 12b or Step 12d does not affect 

50 the result of calculation. In addition, the phase 8 
corresponds to a phase of the generator rotor, and 
of course the same effect is obtained even when 
the phase of the generator rotor is corrected. 

Further, it is possible to empirically determine 

55 the value of n in 360n, but when predetermined 
accuracy and a significant number of digits proper 
to the computer are known, the n can be deter- 
mined as follows. 

10 
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For example, assuming that the significant 
number of digits of the computer is a digits in 
decimal notation and phase accuracy is expressed 
by b digits below the decimal point in a unit of 
degree, a number of digits above the decimal point 
becomes "a - b". Therefore, an integer of "a - b" 
digits becomes the maximum numerical value with- 
in this range. The "a - b" is becomes a number 
smaller than 10 to the "a - b" power. That is, 360n 
should be selected so as to become smaller than 
this value. In the above explanation, the unit Is a 
degree but when using a radian as the unit, the 
predetermined value should be accordingly 
changed to 27rn, and the number of digits b below 
the decimal point should be increased. 

In the same way as the conventional simulator, 
this embodiment is provided with means, that is, 
phase corrector 12 which decides whether the 
phase 5 exceeds a predetermined set value or not 
after calculation of the phase '5 , and which, when It 
exceeds the set value, so controls the phase 5 that 
it always remains within the range of the set values 
by adding or subtracting the predetermined set 
value to or from the phase 6 . Therefore, such 
condition is avoided that the phase 5 falls within an 
error range of the computer system numerical cal- 
culation and thus continuation of calculations be- 
comes impossible. Therefore, the simulation does 
not become Impossible- 
Fig. 1 1 is a block diagram showing a configura- 
tion of another embodiment of the third invention. 
This embodiment includes a plurality of generators 
in the electric power system. In Hg. 11. reference 
numerals 1a. lb ... 1m denote dynamic characteris- 
tic calculators for the respective generators #1. #2 
... #m. Each of these generator dynamic char- 
acteristic calculator la, lb ... 1m has the same 
function as the phase corrector 12 of the above 
embodiment. 

Calculation procedures are as follows. 
First, output voltage Pgii Pg2 Pom and termi- 
nal voltages Vn, V,2 ... V tm of the generators are 
calculated by electric power system network cal- 
culator 2, then the calculated values are inputted to 
the generator dynamic characteristic calculators la, 
lb ... 1pi to calculate amplitudes voi, V02 ... V on, 
and phases 5i, S2 ... 5m of the internal induced 
voltages of the generator. At this time the calcula- 
tions of the generator dynamic characteristic cal- 
culators la, lb ... 1m are not restricted by order of 
the calculations. In the case where there are a 
plurality of calculation units, the calculations may 
be parallel-processed by these units or may be 
serial-processed by one unit 

Fig. 12(a), (b) and (c) are waveform charts 
showing changes in time series of the phases 5i , 
S2 ... 6m of the internal induced voltages of the 
generator, where the abscissas represent times 
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and the ordi nates represent phases (degree) re- 
spectively. 

The phase 6 i shown in Fig. 12(a) tends to 
increase monotonically and when it reaches a pre- 

5 determined values of 360n, it is subtracted by 360n 
to become 0. The phase S 2 shown in Fig. 12(b) 
somewhat changes but does not exceed the range 
of two set values and is not corrected. The phase 5 
3 shown in Fig. 12(c) tends to decrease and when 

70 it reaches a predetermined value of -360n, it is 
subtracted by -360n. that is, added by +360n to 
become 0. 

In this way, the phase 5 is kept within the range 
of -360n to +360n. The third invention is effective 

75 to simulation of a system lasting for a long time 
after the phase 5 has been kept within the range of 
two predetermined set values and a generator has 
caused step-out, or to simulation of a system in- 
cluding a stopped generator, or to simulation of a 

20 system including two A.C. systems with different 
frequencies which are connected with a D.C. trans- 
mission system. 

As mentioned above, the third invention com- 
prises the generator dynamic characteristic calcula- 

25 tor 1 as the first calculating means which calculates 
the generator dynamic characteristic and the elec- 
tric power system network calculator 2 as the sec- 
ond calculating means which calculates power 

' - transmission condition of an electric power system 

30 network. Then the generator dynamic characteristic 
calculator 1 receives information on an output pow- 
er and terminal voltage of the generator outputted 
from the electric power system network calculator 
2 to calculate an amplitude and a phase of the 

35 internal induced voltage of the generator after ev- 
ery one time step of integration to output them to 
the electric power system network calculator 2. and 
the electric power system network calculator re- 
ceives the amplitude and the phase of the internal 

40 induced voltage of the generator outputted from the 
generator dynamic characteristic calculator and 
various constants of the electric power system net- 
work to calculate the power transmission condition 
of the electric power system network at a predeter- 

45 mined time to output information on the output 
power and terminal voltage of the generator to the 
generator dynamic characteristic calculator 1. Fur- 
ther when a phase in the generator dynamic char- 
acteristic calculator exceeds a predetermined set 

50 value obtained by multiplying 360° by an integer, 
a value which is subtracted from this phase value 
by the predetermined value is determined as the 
phase, and then by making the generator dynamic 
characteristic calculator and the electric power sys- 

55 tern network calculator perform calculations alter- 
nately, dynamic changes of the electric power sys- 
tem including generators is serially simulated. 
Thus, in the third invention, after the phase 5 is 
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calculated in the same way as the conventional, a 
check is nnade whether the calculated phase value 
exceeds predetermined set values or not. When it 
exceeds the set value, adding or subtracting the 
set value to or from the phase 5 can always bring 
back the phase 5 within the range of two set 
values, thereby causing no trouble with calculation 
in the computer system. 

Now, the fourth invention will be explained. The 
simulator of the fourth invention comprises the 
same configuration as the simulator which uses the 
conventional decentralized computer shown in Fig. 
2. and calculation procedure itself is basically the 
same as that of the conventional electric power 
system simulator shown in Fig, 3. 

The electric power system simulator shown in 
Fig. 2 uses the decentralized type computer having 
a plurality of CPUs 61, to each of which each 
different equation representing a generator, excita- 
tion system, PSS. or governing system each of 
which is a element of the electric power system is 
assigned. Calculations on the respective elements 
are decentralized-processed by the respective 
CPUs 61, thereby a load on each CPU 61 is 
alleviated and the calculations are performed at a 
high speed as a whole. 

By the way, reference numeral 62 represents a 
communication line " connecting between the re- 
spective CPUs. ^ 

The electric power system, simulator of the 
fourth invention obtains a solution from the follow- 
ing equations (4), (5) and (6) in the network calcula- 
tion shown in Fig. 3. 

Vii(K.i) ^ <^lYo + * Vi (K))A1i (4) 

lii{K + 1) = Yg * Vg + EYij * Vj (k-d 
+ Yii " Vii(K+i) (5) 



liii 



-* 0 



where 

1: current 

Y : admittance 

V: voltage 

That is, the solution employs a method of 
substituting V as a known quantity in an equation 

1 = YV 

and of calculating mere algebraic equation. First, 
using a determinant of 

I = YV. 

The equation (4) is calculated to obtain a new 
Vjifk^-D'The underliried part Yq * Vq of the equation 



(4) is a term corresponding to an equivalent current 
source in the conventional method and is a current 
in a generator bus. In addition, EYij * Vi ^k) is a sum 
of currents in a bus j directly connected to a bus i 

5 at a present time (k). 

Next, using Vij(k*i) representative of new (k+1) 
condition obtained from the equation (4), the equa- 
tion (5) is solved to obtain new lii(k+i)- The equation 
(6) shows that one value of Iii(k + i) obtained from 

70 the equation (5) does not coincide with the other, 
thus the focusing calculation is performed until I 
becomes 0. 

As mentioned above, in the fourth invention, 
when performing calculation concerning one bus. 
75 the calculation is executed using only information 
concerning buses directly connected to the bus 
through a transmission line or a transformer, thus 
the CPUs which simulate the buses should be 
connected each other according to system con- 
20 struction,' as a result focusing of the network cal- 
culation can be carried out at a high speed. Fur- 
ther, provision of a slack bus used for the conven- 
tional bus is not necessary, and because all the 
buses are equal, a bus or transmission line can be 
25 easily added or changed, thus flexibility and 
trouble-withstandability of the decentralized-CPU 
computer system are effectively utilized. 

-In addition.-when-employing a method of cal- 
culation- by using already-obtained V]{k*i) serially 
30 instead of Vj(ic) of equation (4) and obtaining Vi(k^ ^ 
of remaining buses, a number of times of calcula- 
tion until focusing is further reduced. 

As mentioned above, according to the fourth 
invention, when simulation of an electric power 
35 system is performed by the decentralized type 
computer, the network calculation of one bus is 
executed on the basis of only information on buses 
which are directly connected to the bus. Thereby 
the network calculation is fast focused and thus 
40 high speed calculation is enabled and a high speed 
simulator is realized which can analyze an electric 
power system network phenomenon at almost the 
same speed as the real system phenomenon pro- 
ceeding speed. 
45 The electric power system simulator of the fifth 

invention will be explained with reference to a block 
diagram of Fig. 13. 

in Fig. 13 reference numeral 51 denotes the 
same electric power system simulator as a digital 
50 type electric power system simulator utilizing the 
conventional decentralized computer system shown 
in Fig. 2. The fifth invention comprises the digital 
type electric power system simulator 51, an analog 
type electric power system simulator 53 and an 
55 A/D-D/A converter 56 which transfers data between 
both the simulators. 

The analog type electric power system simula- 
tor 53 comprises a plurality of combined miniature 
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simulators 54 each of which simulates each ele- 
ment such as a generator, control system, trans- 
former, load or a relay of an electric power system. 

The digital type electric power system simula- 
tor 51 performs the same calculation as the electric 
power system sinnulator utilizing, for example, the 
decentralized type computer system of the above 
fourth invention. 

Such a digital type electric power system sim- 
ulator 51 and an analog type electric power system 
simulator 53 are connected with the A/D-D/A con- 
verter 56, whereby addition or a change of a sys- 
tem to be simulated is facilitated, as a result a 
detailed analysis of a phenomenon generated in 
the analog type electric power system simulator 53 
becomes possible. 

That is, analog data outputted from the analog 
type electric power system simulator 53 is con- 
verted into digital signals at the A/D-D/A converter 
56 to give them to the computer, whereby phenom- 
ena of each system in the digital type electric 
power system simulator 51 corresponding to each 
simulator 54 configurating the analog type electric 
power system simulator 53 become analyzable. 
Thereby modeling of each system in the digital 
type electric power system simulator 51 can be 
performed more precisely. 

On the contrary, digital data outputted from the 
digital type, electric power system simulator 51 is 
converted into analog signals at the A/D-D/A con- 
verter 56 to add them to analog data outputted 
from the analog type electric power system simula- 
tor 53, whereby an analysis of a large scale electric 
power system which has been impossible with the 
conventional simulator becomes possible with sim- 
ple condition setting. 

Further, analog data outputted from the analog 
type electric power system simulator 53 is substi- 
tuted by analog data outputted from the analog 
type electric power system simulator 53, whereby 
setting of the electric power system can be easily 
changed, in addition initial condition setting of the 
analog type electric power system simulator 53 is 
given from the digital type electric power system 
simulator 51, whereby a desired test condition can 
also be quickly set. 

As mentioned above in detail, according to the 
fifth invention, the analog type electric power sys- 
tem simulator 53 and the digital type electric power 
system simulator 51 which can perform real-time 
processing are connected with the A/D-D/A con- 
verter 56, thus phenomena of a large scale electric 
power system can be simulated by the digital type 
electric power system simulator 51 with simple 
setting. Further not only setting of the electric pow- 
er system is easily changed but also phenomena 
of real machines and apparatuses which are dif- 
ficult to analyze by a computer becomes analyz- 



able with more preciseness, thus simulation of the 
large scale electric power system making the best 
use of both the advantages of the analog type 
electric power system simulator and the digital 
5 type electric power system simulator becomes fea- 
sible with more preciseness and in real time. 

Fig. 14 is a schematic diagram showing a 
configuration of the sixth invention of the electric 
power system simulator, and Fig. 15 is a functional 

70 block diagram of the simulator shown in Fig. 14. 

In the sixth invention, a man-machine device 
57 is connected to the electric power system sim- 
ulator 51 which is the same as the digital type 
electric power system simulator utilizing the con- 

75 ventional decentralized type computer system 
shown in the Fig. 2. 

In the sixth invention as aforementioned, re- 
sults of decentraiized-processed by the digital type 
electric power system simulator 51 are received by 

20 the man-machine device 57 and outputted there- 
from with a technique which depends on human 
sensitivity such as meter indication or graphic dis- 
play. Thus, by watching a man-machine device 57 
screen, an operator can recognize a trouble as if it 

25 were arising just now with meter pointers deflecting 
and trouble data graphic displaying. Further, a con- 
dition change of each apparatus as a corrective 
action to a trouble during simulation can be di- 
rected from the man-machine device 57 as shown. 

30 with an arrow B of Fig. 15, thus simulation of a 
trouble of the real electric power system can be 
performed with a real sense. 

As mentioned above in detail, because the 
system of the sixth invention is provided with the 

35 man-machine device 57. interruption for a condition 
change is enabled at any timing, and condition of 
the electric power system during simulation is dis- 
played with a real sense, and further a setting 
change of condition can be performed during sim- 

40 ulation, thus operation like doing real troubleshoot- 
ing becomes possible. 

As this invention may be embodied in several 
forms without departing from the spirit of essential 
characteristics thereof, the present embodiment is 

45 therefore illustrative and not restrictive, since the 
scope of the invention is defined by the appended 
claims rather than by the description preceding 
them, and all changes that fall within the metes and 
bounds of the claims, or equivalence of such metes 

50 and bounds thereof are therefore intended to be 
embraced by the claims. 

Claims 

55 1. An electric power system simulator compris- 
ing: 

first calculating means for calculating a 
dynamic characteristic of a generator by 
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cyclically obtaining an amplitude and a phase 
of internal induced voltage of said generator 
based on an output power and a terminal volt- 
age of said generator; 

second calculating means for calculating a 5 
dynamic characteristic of a load by obtaining 
an equivalent admittance of said load after 
simulation of the electric power system based 
on an amplitude of an A.C. voltage at an in- 
stallation point of said load; and w 

third calculating means for calculating 
power transmission condition of the electric 
power system network by obtaining the output 
power and the terminal voltage of said gener- 
ator by calculating an A.C. current flowing in 75 
each element of the electric power system with 
using the amplitude and the phase of said 
internal induced voltage of said generator ob- 
tained by said first calculating means, the 
equivalent admittance of said load obtained by 20 
said second calculating means and an imped- 
ance of the electric power system network; 

whereby a condition change of the electric 
power system including said generator is seri- 
ally simulated by alternately performing the 25 
calculations of said first and said second cal- 
culating means and the calculation of said third 
calculating means. . — — 

2. . The electric power system simulator as. set 30/ 

forth in Claim 1, wherein the calculation of A.C. 
current by said third calculating means is per- 
formed by a system network calculation. 

3. An electric power system simulator compris- as 
ing: 

first calculating means for calculating a 
dynamic characteristic of a generator; 

second calculating means for cailculating a 
power transmission condition of an electric 40 
power system network and for outputting an 
output power and a terminal voltage of said 
generator; and 

communicating means for communicating 
between said first calculating means and said 45 
second calculating means by transferring in- 
formation; 

whereby said first calculating means 
cyclically calculates information on an internal 
induced voltage of said generator representa- so 
tive of condition of said generator based on the 
output power and the terminal voltage of said 
generator outputted from said second calculat- 
ing means, and outputs said information to 
said second calculating means through said 55 
communicating means, 

said second calculating means calculates 
power transmission condition of said electric 



power system network at predetermined time 
with using said information and various con- 
stants of said electric power system network, 
and 

a condition change of the electric power 
system including said generator is serially sim- 
ulated by alternately executing the calculation 
of said first calculating means and the calcula- 
tion of said second calculating means, 

4. An electric power system simulator as set forth 
in Claim 3, wherein said second calculating 
means comprises a calculating unit which cal- 
culates the power transmission condition of 
said electric power system network, and a time 
management unit which makes said calculating 
unit calculate by indicating every cycle for 
calculation. 

5, An electric power system simulator compris- 
ing: 

first calculating means, having 
phase calculating means for calculating a 
phase of internal induced voltage of a gen- 
erator, 

phase correcting means for correcting cal- 
culated phase value by said phase calculat- 
_ . - ing means by- subtracting- from the- calcu- - 
lated phase value by value of 360* multi- 
. plied by an integer or adding : to the cal- 
culated phase value by value of 360* multi- 
plied by an integer when an absolute value 
of the phase value calculated by said phase 
calculating means exceeds 360 , and 
amplitude calculating means for calculating 
an amplitude of internal induced voltage of 
said generator, and for calculating a dy- 
namic characteristic of a generator; and 
second calculating means for calculating a 
power transmission condition of an electric 
power system network; 

whereby said first calculating means 
cyclically calculates an amplitude and the 
phase of internal induced voltage of said gen- 
erator representative of condition of said gen- 
erator based on the output power and the 
terminal voltage of said generator outputted 
from said second calculating means, and out- 
puts an information of calculated result to said 
second calculating means, 

said second calculating means calculates 
power transmission condition of said electric 
power system network at predetermined time 
with using said information and various con- 
stants of said electric power system network, 
and 

a condition change of the electric power 
system including said generator is serially sim- 
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ulated by alternately executing the calculation 
of said first calculating means and the calcula- 
tion of said second calculating means. 

6. An electric power system simulator comprising 5 
a decentralized computer including a plurality 
of CPUs and communication lines connecting 
among said plurality of CPUs, by assigning 
each differential equation representative of 
each element constituting electric power sys- w 
tern to each of said plurality of CPUs, which 
simulates phenomena of the electric power 
system with transferring each calculated result 
of sard differential equation decentralized-pro- 
cessed by the respective CPU via said com- is 
munication line, wherein 

the networl< calculation is executed using 
only information on the bus directly connected 
to respective buses when performing the cal- 
culation of a bus being one of said elements. 20 

7- An electric power system simulator compris- 
ing: 

an analog type simulator which simulates a 
phenomenon of an electric power system by 25 
combining a plurality of simulation devices 
each of which simulates each element con- 
stituting electric power system; 

a digital type simulator which calculates a 
phenomenon of the electric power system at a 30 
speed same as or near the proceeding speed 
of the real phenomenon and simulates the 
phenomenon in real time; and 

an A/D-D/A converter which transfers data 
between said analog type simulator and said 35 
digital type simulator. 

8. An electric power system simulator compris- 
ing: 

a digital type simulator which calculates a 40 
phenomenon of an electric power system at a 
speed same as or near the proceeding speed 
of the real phenomenon and simulates the 
phenomenon in real time; and 

a man-machine device which displays in- 45 
formation outpiiited from said digital type sim- 
ulator and is able to perform interruption for a 
condition change to said digital type simulator 
at any time. 
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@ Electric power system simulator. 

@ An electric power system simulator comprising 
first calculator which calculates a generator dynamic 
characteristic, second calculator which calculates a 
load dynamic characteristic and third calculator 
which calculates power transmission condition of an 
electric power system network, then the first calcula- 
tor cyclically obtains an amplitude and a phase of an 
internal induced voltage of the generator based on 
data on an output power and a terminal voltage of 
the generator, and the second calculator calculates a 



load equivalent admittance at predetermined time 
based on data on an A.C. voltage at a load installa- 
tion point, and the third calculator calculates a cur- 
rent flowing in each element of the electric power 
system based on respective data obtained by the 
first calculator and the second calculator and an 
electric power system network impedance to obtain 
the output power and terminal voltage of the gener- 
ator, so that a condition change of the electric power 
system including the generators is serially simulated 
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third calculator. 
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